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ABSTRACT 


This  thesis  is  an  analysis  of  the  one-on-one  ASW  search 
problem  using  a  random  active  search  strategy  in  an  environ¬ 
ment  that  favors  the  target's  counterdetection  ability.  The 
objective  is  to  determine  an  optimum  ping  strategy  by  simu¬ 
lation  of  the  definite-range  problem,  approximation  by  an 
analytical  model  and  use  of  empirical  regression  techniques. 
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INTRODUCTION 


A.  THE  THESIS  OBJECTIVES 

This  thesis  documents  the  analysis  of  one-on-one  ASW 
encounters  between  a  surface  searcher  using  active  sonar  and 
an  evasive  target  submarine.  The  analysis  is  based  on  data 
generated  by  a  computer  simulation  of  the  relative  motion  of 
the  two  adversaries  over  time.  The  specific  objective  of 
this  analysis  is  to  prescribe  a  strategy  for  selecting  a 
searcher  ping  interval  which  maximizes  the  probability  of 
detection  in  an  environment  which  favors  the  target's  counter 
detection  ability. 


B.  THE  SCENARIO 

A  single  surface  shi.p  is  assigned  to  search  for  a  sub¬ 
marine  target  of  interest  using  active  sonar  within  a  region 
several  hundred  thousand  square  miles  in  area.  The  acoustic 
environment  is  considered  homogeneous  throughout  the  area. 
Thus  for  any  particular  case,  the  sonar  range  is  considered 
a  constant.  -  ^ 

1 .  The  Searcher's  Tactic 

The  searcher's  tactic  is  to  move  through  the  area  at 
a  set  speed,  changing  course  randomly  at  times  described  by 
an  exponential  distribution  of  mean  1/0.  The  use  of  the 
exponential  distribution  for  this  purpose  seems  tactically 
prudent  because  of  its  memoryless  property.  The  searcher 


pings  at  times  selected  at  random  from  some  distribution. 

By  this  tactic,  the  exact  time  between  successive  pings  is 
made  unpredictable. 

2 .  The  Target’s  Tactic 

The  target  is  patrolling  the  area  of  interest  at  a 
set  speed,  changing  course  randomly  at  times  also  selected 
from  an  exponential  distribution,  not  necessarily  the  same 
as,  but  independent  of,  that  of  the  searcher.  The  target  is 
capable  of  counterdetecting  the  searcher's  transmissions  at 
ranges  greater  than  the  searcher's  detection  ranges.  It  is 
assumed  that  the  target  has  no  method  of  detecting  and 
locating  the  searcher  other  than  by  passive  detection  of  the 
searcher's  transmissions.  If  the  target  does  counterdetect 
the  searcher  outside  the  searcher's  detection  range,  he 
sprints  away  from  the  searcher  radially  at  a  speed  greater 
than  that  of  the  searcher.  The  duration  of  this  sprint  is  a 
tactical  decision  made  by  the  target,  based  on  what  he  con¬ 
siders  a  "safe"  range,  and  can  only  be  estimated  by  the 
searcher . 
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II .  MODEL  DESCRIPTION 

It  is  commonly  assumed  that  the  time,  Tq ,  required  for 
a  randomly  moving  searcher  and  target  to  first  come  within 
some  relatively  small  range  of  one  another  is  distributed 
exponentially  with  a  mean  that  is  a  function  of  that  range, 
the  searcher's  and  target's  speeds,  and  the  size  of  the  area 
in  which  they  are  confined  [Ref.  1],  Intuitively,  it  would 
seem  if  the  searcher  selects  a  ping  strategy  that  maximizes 
the  probability  of  detection,  given  that  the  target  is 
within  counterdetection  range,  that  strategy  tends  to  mini¬ 
mize  the  expected  total  time,  T,  spent  searching  in  the  area 
for  the  target.  With  this  in  mind,  a  model  of  the  search 
problem  after  the  realization  of  TQ  is  described  below. 

A.  THE  EVENT  DISK 

In  the  model  used  for  simulation,  the  event  disk,  i.e., 
the  region  within  an  event  circle  of  radius  C,  the  counter¬ 
detection  range,  represents  the  area  in  which  any  interaction 
between  the  target  and  searcher  must  occur.  Concentric  with 
the  event  disk  is  the  searcher's  detection  disk  of  radius  D. 
The  event  disk  is  centered  on  that  opponent  with  the  higher 
speed.  Figure  2.1  illustrates  the  case  for  which  the  searcher 
is  at  a  higher  speed  than  the  target,  but  the  labeling  is 
completely  arbitrary  because  of  symmetry.  If  the  target  is  at 


■*.  >>» 


the  higher  speed,  he  is  placed  in  the  center.  All  relative 
relationships  remain  the  same. 


B.  THE  EXPECTED  SEARCH  TIME 

Once  the  target  has  entered  the  event  disk,  one  of  three 
events  must  occur: 

(1)  The  target  departs  the  event  disk  before  the  searcher’s 
first  ping,  by  virtue  of  the  relative  motion  between 
the  two . 

(2)  The  target  sprints  out  of  the  event  disk  as  a  result 
of  being  located  in  the  counterdetection  zone  but  not 
in  the  searcher's  detection  zone  when  the  searcher 
pings . 

(3)  The  target  is  detected  as  a  result  of  being  located  in 
the  detection  zone  when  the  searcher  pings.  This  event 
completes  the  search. 

If  Events  (1)  or  (2)  occur,  then  there  is  a  possibility  that, 
eventually,  the  target  will  again  enter  the  event  disk. 
Therefore,  once  the  target  has  entered  the  event  disk  for  the 
first  time,  the  remainder  of  the  search  can  be  thought  of  as 
a  series  of  cycles  during  which  the  target  is  either  detected 
or  not  detected.  This  suggests  the  use  of  the  geometric  dis¬ 
tribution  to  describe  the  process.  If  the  search  requires  N 
s,.ch  cycles,  then  (N-l)  of  the  cycles  must  have  resulted  in 
no  detection  occurring.  Therefore,  if  P  is  the  probability  of 
detection,  the  probability  that  the  search  requires  n  cycles 
is 


P  (N  =  n ) 


(l-P)n_1P 


(2.1) 


Determining,  E[T],  the  expected  total  time  for  completing 
the  search  requires  that  two  additional  variables  be  defined. 


Let 


Tc (i)  i  =  1 , 2 , 3 ,  .  .  .  , n 

be  a  sequence  of  independent,  identically  distributed  random 
variables  describing  the  cycle  time  ror  the  ifc^  cycle  and 
let  be  the  time  required  for  the  target  to  enter  the  detec¬ 
tion  zone  and  be  detected,  given  he  is  located  on  the  perimeter 
of  the  event  disk.  Then 


N-l 

T  =  Tn  +  l  T  (i)  +  T  (2.2) 

U  i=l  C  a 


and 


E[T|N]  =  E[Tq]  +  (n+l)E[TcJ  +  E[Td]  (2.3) 


Removing  the  condition  on  n  results  in  the  following 
expression : 


co 

E  [T]  =  E  { E  [T  [  N)  }  =  E[Tq]  +  E  [Tj  £  (n-1 )  ( 1-P) n_1P 

n=l 


+  E[TdJ 


(2.4) 


10 


i  —  »  .J 


The  summation  term  is  easily  collapsed.  Let 


S  =  l  (n-1 ) (1-P) n  1P  =  l  n (1-P) nP 
n=l  n=0 

=  0  +  (l-P)P  +  2 (1-P) 2P  +  ...  (2.5) 

Then 

(1-P) S  =  (1-P)2P  +  2(1-P)3P  +  3  ( 1-P) 4P  +  ...  (2.6) 

Subtracting  Equation  (2.6)  from  Equation  (2.5)  yields 

PS  =  P(l+P)[l  +  (1+P)  +  (1+P)2  +  ...]  (2.7) 

The  sum  inside  the  brackets  is  a  geometric  series  which  con¬ 
verges  to  1/P.  Therefore 

E  [T]  =  E(TQ]  +  (i-l)E[Tc]  +  E[TdJ  (2.8) 

It  can  be  seen  that  maximizing  the  probability  of  detection, 
P,  will  aid  in  minimizing  the  expected  total  search  time.  It 
is  for  this  reason  that  this  thesis  concentrates  on  the  prob¬ 
lem  within  the  event  disk;  that  is,  finding  a  ping  rate  which 
maximizes  P. 
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C.  CRITICAL  ASSUMPTIONS 


Several  assumptions  have  been  made  to  simplify  the  model 
and  analysis  of  the  generated  data.  In  addition  to  those  stated 
previously,  the  following  also  apply: 

(1)  The  occurrence  of  detection  and  counterdetection  events 
is  determined  using  a  definite-range  or  "cookie 
cutter"  model. 

(2)  The  target  and  searcher  have  negligible  baffle  areas. 

(3)  Detection  and  counterdetection  ranges  are  not  degraded 
with  speed. 

(4)  There  is  no  convergence  zone  considered,  nor  are  there 
any  gaps  in  the  event  circle. 

(5)  The  target  is  strictly  evasive. 

(6)  Counterdetection  range  as  used  in  the  model  should  be 
considered  t'..a  target's  minimum  desired  range  to  the 
searcher.  It  is  assumed  this  range  is  at  least  twice 
the  detection  range. 


D.  THE  REQUIREMENT  FOR  A  PING  STRATEGY 

The  assumed  existence  of  a  definite  counterdetection  range 
greater  than  the  searcher's  detection  range  requires  that  the 
searcher  have  a  well-defined  minimum  interval  between  any  two 
active  pings.  This  minimum  interval  is: 


I 


min 


C  -D 
U  +  V 


(2.9) 


where  C  is  the  counterdetection  range,  D  the  detection  range, 
and  the  denominator  is  the  sum  of  the  two  speeds.  This  is 
merely  the  time  required  for  the  target  to  move  from  the 
perimeter  of  the  event  disk  to  the  detection  zone  at  the  maxi¬ 
mum  attainable  relative  speed. 


There  is  also  a  maximum  practical  ping  interval  that  is 
not  as  well  defined.  That  is,  if  the  searcher  pings  very 
infrequently  he  loses  the  opportunity  for  detection  because 
the  target  may  transit  in  and  then  out  of  the  detection  zone 
between  pings.  This  "maximum"  should  depend  upon  the  size  of 
the  detection  zone,  and  the  relative  motion  between  the 
searcher  and  target.  It  is  because  of  the  relative  motion 
aspects  that  this  maximum  practical  interval  cannot  be  defined 
as  easily  as  the  minimum.  The  existence  of  a  minimum  ping 
interval  below  which  the  probability  of  detection  is  zero, 
and  a  "maximum"  ping  interval  beyond  which  the  probability  of 
detection  is  small,  implies  that  the  probability  of  detection 
reaches  a  maximum  between  these  two  extremes.  This  maximum 
should  be  a  function  of  the  ranges  and  speeds  specific  to  each 
particular  case.  Therefore  the  first  step  is  to  determine  a 
ping  strategy  as  it  depends  upon  the  independent  variables. 

This  will  be  done  by  simulation,  approximate  analytical  modeling, 
and  a  blending  of  the  two  by  an  empirical  regression  technique. 

E.  DATA  GENERATION 

The  equation  for  probability  of  detection  if  both  searcher 
and  target  remain  on  constant  courses  and  speeds  is  complex 
but  can  be  solved  using  polar  coordinates  and  some  trigonometry. 
However,  in  the  problem  at  hand,  both  relative  speed  and  its 
direction  change  randomly.  This  urges  the  use  of  simulation 
to  generate  data  on  relative  courses  and  positions. 


The  simulation  program  for  the  definite  range  problem 
(included  in  Appendix  C  for  informational  purposes)  supplies 
as  output  the  number  of  detections,  the  number  of  counterde¬ 
tections,  and  the  number  of  times  the  target  departs  the  event 
disk  before  the  first  ping,  for  a  given  ping  interval,  detec¬ 
tion  range,  counterdetection  range,  searcher  speed,  and  target 
speed.  The  frequency  of  course  changes,  the  searcher  course, 
and  the  target  course  are  determined  by  random  number  generators 
All  relative  motion  is  placed  on  the  target  using  trigonometric 
arguments  with  the  searcher  remaining  at  the  coordinate  origin. 

Because  each  interaction  begins  with  the  range  between 
the  searcher  and  target  decreasing  to  C,  it  follows  that  the 
initial  relative  velocity  must  be  directed  into  the  event 
disk.  To  accomplish  this  in  the  simulation,  the  first  rela¬ 
tive  velocity  vector  was  determined  by  assuming  that  the 
searcher's  speed  component  directly  toward  the  target  was  just 
greater  than  the  target's  speed.  After  the  initial  leg,  all 
motion  was  unconstrained.  So  after  several  course  changes, 
the  effect  of  the  initial  leg  is  "forgotten"  by  the  process. 

An  alternate  method,  not  used  in  this  thesis,  would  be  to 
let  the  target's  initial  relative  angle  on  the  bow  be  selected 
from  a  cosine  distribution.  That  is, 

l  4> 

P  (Relative  AOB  <_<}>)  =  y  /  cos  9  d8 

1  -tt/2 


i(sin  4)  +1),  -tt/2  <  <p  <_  tt/2  (2.10) 


This  procedure  is  suggested  by  Koopman's  observation  that 
when  searching  for  a  stationary  target,  the  bearing  of 
initial  detection  will  have  a  cosine  distribution  [Ref.  1] . 

It  is  not  known  which  of  these  methods  is  more  correct. 

It  may  be  possible  to  derive  an  exact  expression  for  the  joint 
distribution  of  the  target's  relative  speed  and  course  at 
detection,  but  this  was  not  accomplished  here. 

Once  the  time  for  the  first  ping  is  reached,  the  replica¬ 
tion  is  stopped.  The  range  between  the  two  is  calculated 
and  compared  to  the  values  for  C  and  D  to  determine  which  of 
the  possible  events  has  occurred.  The  outcome  is  stored  and 
the  whole  process  is  repeated  for  the  desired  number  of  trials 
The  data  is  then  analyzed  to  obtain  estimates  of  the  proba¬ 
bility  of  detection  and  other  relevant  quantities,  such  as 
expected  times  within  and  without  the  event  disk. 


III.  DETERMINATION  OF  THE  REGRESSION  MODEL 


Once  the  data  has  been  generated,  regression  may  be  used 
to  determine  an  empirical  predictive  formula  for  the  optimum 
ping  rate.  Armed  with  only  intuitive  hypotheses,  the  search 
for  the  "best"  functional  form  of  the  input  variables  would 
be  difficult.  So  it  is  desirable  to  find  some  theoretical 
guidelines  for  selecting  candidate  explanatory  variables  to 
use  in  the  regression  model. 

A.  USE  OF  THE  VON  NEUMANN  FUNCTION 

The  theoretical  model  selected  for  use  is  one  for  energy 
transmission  and  return  when  the  target's  motion  is  modeled 
by  a  diffusion  process.  Define  the  searcher's  location  as 
the  origin  on  a  Cartesian  plane,  and  define  the  target's 
location  at  time  t,  as  [X(t),Y(t)].  Let  the  target  undergo 
Brownian  motion  so  that  its  location  at  time  t  is  described 
by 


X(t)  ~  N  (X  (0  ,cr2t) 

Y  (t )  ~  N  (Y  ( 0 )  , a  2t ) 


(3.1) 


and  let  the  searcher's  probability  of  detection,  for  a  ping  at 
time  t  be 
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P (t | X ( t) , Y (t ) )  =  e-(X(t)2+Y(t)2)/«2 


=  e 


1  2  2 
-j(R(tr)/6z 


(3.2) 


where  R(t)  is  the  range  from  the  searcher  to  the  target  and 
.2 


<S  is  a,  thus  far,  unspecified  constant.  The  constant  6  in 
the  detection  function  (Equation  (3.2))  plays  a  role  analogous 
to  that  of  the  detection  disk  radius,  D.  In  particular,  6  is 


-1/2 

that  range  where  the  probability  of  detection  is  e  ~  0.607, 


Removing  the  condition  on  position  by  integrating  over 
all  values  of  X(t)  and  Y(t)  results  in  the  following  expression 
for  P  (t)  : 


-  -y (X (t ) 2 ) /S  2  .  -it ( (X(t)-X(O) )2)/o2t] 

P(t)  =  /  /e  - - -  e  ^  dX 


'J  2  tt  o  2 1 


-j(  (Y(t)2)/62  -j[ ( (Y(t)-Y(O) )2)/a2t] 

e  -  e  dY 


v/27T02t 


=  e 


-j[(R(0)2/62)/(l+o2t/fi2) ] 


2  2 
(i  +o  t/6  : 


(3.3) 


Using  the  natural  log  function  to  linearize  P(t),  and 
differentiating  with  respect  to  t  yields 


d  HnP(t) 


2  /x2 
o  /6 


1  R 


dt 


[±.  ~ _ 

l  o  O 


~2  5" 1 2  2 

1  +0  t/6  6 i 


(1  +q~^)  1 


(3.4) 
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Setting  the  derivative  equal  to  0  results  in  the  following 
expression  for  T*,  the  optimum  ping  interval: 


T* 


1,R(0)2' 

*  r 


1 


Therefore 


T* 


0  if 


1  R ( 0)  2  ^ 

2  62  " 


1 


and 


T* 


1,R(0) '2 
2k  6  ' 

- J? 

V 


1 


otherwise 


(3.5) 


(3.6)  ; 


Referring  to  Equation  (3.2),  suppose  energy  from  a  ping 
transmitted  from  location  (0,0)  reaches  the  target  with  ( 

4 

probability  . 


1  ,R ( t) , 2 

2  6o 


I 


(3.7) 


and  the  probability  of  the  signal  returning  to  the  searcher 
given  it  reached  the  target  is 


1 <R ( t) ,2 
e  '' 
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Then  the  probability  of  the  searcher  detecting  the  target 
becomes 


_1/R(t),2  _l,R(t),2 

2"  ~ 5 - '  2"  <5~ 

P (detection  R(t))  =  e  °  e  1 

-i-R  ( t )  2  [~\  +  4*-] 

5^  6*: 

o  1 

=  e 

lfR(t)  2 

=  e'2  6  (3.9) 

which  is  also  a  Von  Neumann  detection  function. 

The  intention  at  this  point  is  to  use  the  form  of  Equa¬ 
tion  (3.5)  to  develop  an  approximation  for  the  optimal  ping 
interval  T*  for  the  definite  range  problem.  This  is  accom- 
plished  by  setting  R(0)  in  Equation  (3.5)  equal  to  the 
counterdetection  range,  C,  and  recognizing  that  the  parameter 
6  is  a  "characteristic  range"  for  the  Von  Neumann  function. 

So 


1  -  r  1  +  1  1 

7*  ~  [7?  +  72] 
6  6q  6i 


[-T  +  ~T] 

R  r 


(3.10) 


where  R  and  r  are  the  ranges  associated  with  the  outbound 
and  inbound  acoustic  paths. 

Substituting  into  Equation  (3.5)  yields 


T* 

x  no 


\c 2  (“7  +  '  1 

C  r 


(3.11) 


Letting  r  now  be  the  active  detection  range,  D,  results  in 
the  following  expression: 


t4,  tvc.2  -_yp^  2/c2, 

2 a  1/C  +  1/D^ 


2  2  2 

r  ,  fl/D^  -  1/c 
L  2  J  I  o  5" 

2a  1/D^  +  1/C 


] 


C 

2a' 


C/D  +  1 


(3.12) 


B.  THE  EXPLANATORY  VARIABLES 

It  remains  to  replace  the  diffusion  constant  o  with  appro¬ 
priate  random  tour  model  parameters  u,  v,  and  A.  Lambda  is 
the  parameter  for  the  exponential  distribution  describing  the 

minimum  time  between  course  changes  for  either  of  the  adver- 
2  . 

saries.  If  a  is  a  diffusion  rate  describing  the  relative 
motion  of  the  two  then 


a 


2 


U2  + 


(3.13) 


is  dimensionally  correct  and  has  some  theoretical  appeal. 

Specifically,  an  unconstrained  two  dimensional  random  tour 

with  constant  speed  V  and  rate  of  course  change  Av  can  be 

approximated  for  large  times  by  a  diffusion  process  with 
2 

constant  V  /A^lRef.  2].  When  two  particles  are  simultaneously 


■  m  t  ^  CTrwrp  wwwfc  ■!'  p  l  ■  t  ■  im  11  gf  wi  i 


■■WW 


conducting  random  tours,  then  the  composite  random  tour  in 
relative  space  has  a  rate  of  course  change: 


X  =  A  +  X  (3.14) 

u  v 

and  a  random  speed  SR.  If  the  angle  between  the  V  and  U 

velocity  vectors  is  uniformly  distributed  between  0  and  2ir 

,  2  . 

radians  then  the  expected  value  of  SR  is: 


E[S2] 


1  2  7T  n  a 

^  i  (U  +  V  -2  UV  cos  0 )  d0 
2t i  > 


2  2 
U  +  V 


(3.15) 


Using  (U2  +V2)  as  a  representative  squared  speed  for  the 
composite  random  tour  and  X  for  the  rate  of  course  changes 
yields  Equation  (3.13). 

Equation  (3.12)  is  related  to,  but  probably  unequal  to, 
the  optimum  ping  interval  for  the  definite  range  law  model 
that  is  being  simulated.  In  order  to  better  adapt  the  Von 
Neumann  model  results  to  the  definite  range  law  data,  one 
can  redefine  T*  as 


T* 

DR 


31C 


AC' 


(U~  +V“) 


62fC2/D2  -  1j33 

T7]  1  2  '  J 


C  /D  +  1 


(3.16) 


and  determine  the  parameters  8^,  82  >  an<^  S3  by  regression, 
using  as  explanatory  variables,  the  quantities 
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It  is  acknowledged  that  Equation  (3.16)  is  probably  not 
the  "best"  definite  range  law  extrapolation  of  the  Von  Neumann 
result  of  Equation  (3.5).  In  particular,  setting  6^  equal 
to  the  detection  range  D  seems  suspect  since  the  active 
detection  process  involves  two-way  propagation  and  6^  considers 
primarily  the  return  path.  Nonetheless,  Equation  (3.16)  was 
tested  as  a  candidate  regression  model  and,  as  the  next  section 
documents,  it  performed  quite  well. 


IV.  THE  RESULTS 


A  total  of  180  cases  were  run  using  4  different  speed 
combinations,  detection  ranges  from  3  to  16  miles,  and  several 
counterdetection  ranges  between  20  and  40  miles.  The  proba¬ 
bility  of  detection  was  measured  at  each  of  50  different  ping 
intervals,  beginning  at  the  practical  minimum  and  stepped  in 
0.05  hour  increments.  A  total  of  500  trials  were  run  at  each 
ping  interval.  The  observed  T*  was  that  ping  interval  at 
which  the  maximum  probability  of  detection  occurred  for  each 
case.  In  the  event  of  a  tie,  the  earlier  time  was  used. 

Figure  4.1  illustrates  the  dependence  of  the  probabilities  of 
detection,  counterdetection,  and  departure  on  the  value  of 
ping  interval.  As  expected,  the  probability  of  detection  is 
0  until  the  ping  interval  is  greater  than  (C-D)/(U+V),  in¬ 
creases  to  a  maximum  and  slowly  decreases  to  a  small  positive 
value.  The  probability  of  departure  increases  with  ping 
interval,  while  the  probability  of  counterdetection  decreases. 


A.  THE  EMPIRICAL  PREDICTION  FORMULA 

Performing  linear  regression  on  the  data,  using  the 
explanatory  variables  discussed  previously,  produced  the 
following  equation  for  predicting  T*R: 


T* 

DR 


=  0 . 74  [ 


AC 


,0.661 rC2/D2  -1,0.173 

(tr  +v  )  ct/d  +i 


(4.1) 


v.i 


COUNTERDETECTION 


Figure  4.1  Time  Dependency  of  Event  Probabilities 


This  formula  explained  89%  of  the  total  variation  in  the  data 
and  all  coefficients  were  significant  at  the  0.99  level  using 
Student's  t  statistic. 

Tables  ( 1 )  —  ( 4 )  of  Appendix  A  list  the  input  variables, 
observed  T*  and  predicted  T*  for  each  case.  It  is  immediately 
apparent  that  the  probability  of  detection  at  the  predicted 
value  of  T*  is  consistently  less  than  or  equal  to  that  at 
the  observed  value  of  T*.  This  is  because  of  the  initial 
relatively  crude  method  used  to  choose  the  observed  T*  for 
each  case.  The  values  for  the  difference  in  the  probabilities 
a  measure  of  the  prediction  validity,  are  relatively  small. 
Figure  4.2  is  a  scatter  plot  of  the  probability  of  detection 
at  the  predicted  T*  versus  that  at  the  observed  T*.  The  ideal 
•-  lot  would  be  a  straight  line  of  slope  1  through  the  origin. 
The  least  squares  fitted  line  throudi  the  data  has  a  slope  of 
.  }  4  ind  intercept  of  -0.01.  Many  of  the  points  away  from  the 
:  :a  : a  1  can  be  explained  by  examination  of  the  raw  data. 
i ten,  the  predicted  and  observed  optimum  ping  intervals  are 
within  0.10  hours  of  one  another,  but  the  variance  of  the 
sampled  binomial  distribution  causes  the  two  detection  proba¬ 
bilities  to  appear  farther  apart  than  might  actually  be  the 
case.  It  should  be  mentioned  that  smoothing  half  of  the  raw 
data  using  running  medians  followed  by  running  averages 
(Hanning)  and  using  the  same  regression  model  did  not  alter 
the  coefficients  of  the  prediction  formula  significantly. 

It  did,  as  expected,  tend  to  decrease  the  difference  between 


93948 


PROBABILITY  OF  DETECTION  AT  OBSERVED  T* 


the  value  for  probability  of  detection  at  the  predicted  T* 
and  that  at  the  observed  T*. 

B.  CONCLUSIONS 

The  most  significant  finding  is  that  it  is  possible  to 
determine  an  optimum  ping  strategy  within  the  limitations  of 
the  model.  If  a  searcher  were  tasked  with  such  a  search,  he 
would  not  ever  want  to  ping  more  frequently  than  the  minimum 
ping  interval  (C-D)/(U+V).  This  assumes  the  searcher  has  a 
good  estimate  of  the  counterdetection  range  and  speed  of  the 
target.  On  the  average,  he  would  want  to  ping  at  a  rate  that 
is  slightly  less  than  that  prescribed  by  the  prediction. 

There  are  two  reasons  for  this:  (1)  the  probability  of 
detection  decreases  more  slowly  to  the  right  of  T*  than  it 
increases  at  the  left;  and  (2)  any  deviation  to  the  right  of 
T*  merely  increases  the  probability  of  the  departure  event 
occurring.  The  departure  event  is  considered  less  detrimental 
to  the  search  effort  than  the  counterdetection  event,  for 
which  the  probability  of  occurrence  steadily  decreases. 

The  Von  Neumann  function  seems  well-suited  as  a  theoreti¬ 
cal  foundation  for  determination  of  an  optimum  ping  strategy. 
This  suggests  that  research  into  its  use  in  a  more  realistic 
model  of  active  sonar  search  might  prove  valuable  in  predictions 
of  this  sort.  The  extremely  sharp  cut  off  of  the  definite 
range  law  is  certainly  an  artificiality. 
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TABLE  1 

COMPARISON  OF  OBSERVED  AND  PREDICTED  T* 


SEARCHER  SPEED: 

15  KTS 

TARGET 

SPEED: 

5  KTS 

LAMBDA:  1. 

5/HR 

C(NM.  ) 

D(NM.  ) 

OBSERVED 

PREDICTED 

T* 

P(  T* ) 

x* 

P(  T* ) 

20 

3 

1.  75 

o:  030 

1.  30 

0.  022 

4 

1.  30 

0.  058 

1.  30 

0.  058 

5 

1.  55 

0.  070 

1.  30 

0.  040 

6 

1.  35 

0.  094 

1.  30 

0.  056 

7 

1.  45 

0.  146 

1.  25 

0.  126 

8 

1.  25 

0.  156 

1.  25 

0.  156 

22 

4 

1.  65 

0.  042 

1.  50 

0.  018 

5 

1.  90 

0.  056 

1.  45 

0.  040 

6 

1.  20 

0.  076 

1.  45 

0.  066 

7 

1.  25 

0.  122 

1.  45 

0.  068 

8 

1.  40 

0.  130 

1.  45 

0.  124 

9 

1.  75 

0.  162 

1.  40 

0.  158 

25 

4 

1.  80 

0.  028 

1.  75 

0.  020 

5 

1.  35 

0.  058 

1.  75 

0.022 

6 

1.  70 

0.  060 

1.  75 

0.  048 

7 

1.  80 

0.  092 

1.  70 

0.052 

8 

1.  85 

0.  100 

1.  70 

0.096 

9 

1.  85 

0.  128 

1.  70 

0.  112 

10 

1.  75 

0.  140 

1.  70 

0.  134 

27 

4 

2.  05 

0.  028 

1.  95 

0.  006 

5 

2.  60 

0.  040 

1.  95 

0.  028 

6 

2.  00 

0.  054 

1.  95 

0.  028 

7 

2.  35 

0.  064 

1.  90 

0.  058 

8 

2.  00 

0.  084 

1.  90 

0.  078 

9 

1.  55 

0.  098 

1.  90 

0.  086 

10 

2.  10 

0.  132 

1.  90 

0.  076 

11 

1.  95 

0.  136 

1.  85 

0.  106 

30 

6 

2.  90 

0.  042 

2.  25 

0.  026 

7 

2.  10 

0.  060 

2.  20 

0.  030 

8 

1.  90 

0.  066 

2.  20 

0.  050 

9 

2.  40 

0.  078 

2.  20 

0.  070 

10 

2.  00 

0.  092 

2.  15 

0.  074 

11 

2.  35 

0.  146 

2.  15 

0.  106 

12 

1.  85 

0.  126 

2.  15 

0.  112 

13 

2.  20 

0.  162 

2.  10 

0.  112 

14 

2.  00 

0.  164 

2.  10 

0.  142 

40 

8 

3.  70 

0.  044 

3.  25 

0.  034 

9 

3.  30 

0.  052 

3.  25 

0.  042 

10 

3.  10 

0.  046 

3.  25 

0.  042 

11 

2.  80 

0.  064 

3.  20 

0.  038 

12 

3.  65 

0.  066 

3.  20 

0.  046 

13 

3.  05 

0.  080 

3.  20 

0.  062 

14 

3.  05 

0.  082 

3.  15 

0.  068 

15 

2.  85 

0.  112 

3.  15 

0.  076 

16 

2.  35 

0.  112 

3.  10 

0.  082 
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TABLE  2 

COMPARISON  OF  OBSERVED  AND  PREDICTED 

T* 

SEARCHER  SPEED:  18  KTS 

TARGET 

SPEED: 

7  KTS 

LAMBDA: 

1.  5/HR 

C(NM.  ) 

D(NM.  ) 

OBSERVED 

PREDICTED 

T* 

P(  T*  } 

■J* 

P(  T* ) 

20 

3 

1.  05 

0.  026 

1.  00 

0.  020 

4 

0.  90 

0.  044 

1.  00 

0.  034 

5 

1.  00 

0.  074 

1.  00 

0.  074 

6 

1.  00 

0.  092 

1.  00 

0.  092 

7 

1.  00 

0.  130 

0.  95 

0.  122 

8 

0.  95 

0.  190 

0.  95 

0.  190 

22 

4 

1.  30 

0.  040 

1.  15 

0.  028 

5 

1.  25 

0.  054 

1.  15 

0.  040 

6 

1.  15 

0.  078 

1.  10 

0.  060 

7 

1.  35 

0.  098 

1.  10 

0.  076 

8 

1.  00 

0.  122 

1.  10 

0.  106 

9 

1.  25 

0.  178 

1.  10 

0.  132 

25 

4 

1.  30 

0.  032 

1.  35 

0.  016 

5 

1.  40 

0.  054 

1.  35 

0.  022 

6 

1.  50 

0.  062 

1.  35 

0.  044 

7 

1.  40 

0.  086 

1.  30 

0.  076 

8 

1.  60 

0.  096 

1.  30 

0.  084 

9 

1.  55 

0.  138 

1.  30 

0.  110 

10 

1.  25 

0.  168 

1.  30 

0.  148 

27 

4 

1.  30 

0.  024 

1.  50 

0.  012 

5 

1.  40 

0.  040 

1.  50 

0.  032 

6 

2.  00 

0.  060 

1.  50 

0.  034 

7 

1.  20 

0.  080 

1.  45 

0.  050 

8 

1.  70 

0.  096 

1.  45 

0.  090 

9 

1.  50 

0.  112 

1.  45 

0.  090 

10 

1.  75 

0.  132 

1.  45 

0.  094 

11 

1.  35 

0.  160 

1.  40 

0.  144 

30 

6 

1.  80 

0.  040 

1.  70 

0.  016 

7 

2.  00 

0.  056 

1.  70 

0.  038 

8 

1.  80 

0.  078 

1.  70 

0.  056 

9 

1.  70 

0.  082 

1.  70 

0.  082 

10 

1.  60 

0.  110 

1.  65 

0.  088 

11 

2.  15 

0.  128 

1.  65 

0.  108 

12 

1.  70 

0.  154 

1.  65 

0.  148 

13 

1.  65 

0.  168 

1.  60 

0.  130 

14 

1.  50 

0.  196 

1.  60 

0.  174 

40 

8 

2.  50 

0.  042 

2.  50 

0.  042 

9 

2.  65 

0.  048 

2.  50 

0.  022 

10 

2.  45 

0.  058 

2.  50 

0.  050 

11 

2.  35 

0.  068 

2.  45 

0.  042 

12 

2.  20 

0.  078 

2.  45 

0.  052 

13 

2.  45 

0.  088 

2.  45 

0.  088 

14 

2.  15 

0.  100 

2.  40 

0.  084 

15 

3.  25 

0.  108 

2.  40 

0.  074 

16 

2.  20 

0.  120 

2.  40 

0.  090 

TABLE  3 

COMPARISON  OF  OBSERVED  AND  PREDICTED  T* 


C(NM.  ) 
20 


22 


25 


27 


30 


40 


SEARCHER  SPEED:  20  KTS 
TARGET  SPEED:  10  KTS 
LAMBDA:  1. 5/HR 


D( NM. )  OBSERVED  PREDICTED 

'  m  /  m  j.  \  mu.  /  mu. 


1  t 

3 

T* 

0.  75 

P(TM 

0.  036 

T* 

0.  85 

P(TM 
0.  022 

4 

1.  00 

0.  074 

0.  80 

0.  044 

5 

0.  80 

0.  080 

0.  80 

0.  080 

6 

1.  00 

0.  144 

0.  80 

0.  088 

7 

0.  80 

0.  154 

0.  80 

0.  154 

8 

0.  85 

0.  170 

0.  80 

0.  166 

4 

0.  80 

0.  048 

0.  95 

0.  038 

5 

0.  85 

0.  080 

0.  95 

0.  058 

6 

0.  85 

0.  082 

0.  90 

0.  082 

7 

0.  95 

0.  118 

0.  90 

0.  090 

8 

1.  00 

0.  144 

0.  90 

0.  132 

9 

0.  85 

0.  164 

0.  90 

0.  164 

4 

0.  95 

0.  032 

1.  10 

0.  016 

5 

0.  95 

0.  058 

1.  10 

0.  032 

6 

1.  15 

0.  064 

1.  10 

0.  044 

7 

1.  30 

0.  086 

1.  10 

0.  058 

8 

0.  95 

0.  124 

1.  10 

0.  078 

9 

0.  90 

0.  128 

1.  05 

0.  090 

10 

0.  90 

0.  158 

1.  05 

0.  132 

4 

1.  10 

0.  026 

1.25 

0.  014 

5 

1.  10 

0.  038 

1.  20 

0.  024 

6 

1.  00 

0.  062 

1.  20 

0.  042 

7 

1.  10 

0.  088 

1.  20 

0.  040 

8 

1.  10 

0.  086 

1.  20 

0.  074 

9 

1.  25 

0.  120 

1.  20 

0.  076 

10 

1.  40 

0.  126 

1.  20 

0.  104 

11 

1.  20 

0.  162 

1.  15 

0.  142 

6 

1.  20 

0.  040 

1.  40 

0.  028 

7 

1.  80 

0.  062 

1.  40 

0.  034 

8 

1.  30 

0.  078 

1.  40 

0.  050 

9 

1.  35 

0.  090 

1.  40 

0.  068 

10 

1.  45 

0.  108 

1.  35 

0.  084 

11 

1.  50 

0.  132 

1.  35 

0.  102 

12 

1.  30 

0.  164 

1.  35 

0.  148 

13 

1.  25 

0.  176 

1.  35 

0.  158 

14 

1.  25 

0.  188 

1.  30 

0.  158 

8 

1.  55 

0.  038 

2.  05 

0.  020 

9 

1.  70 

0.  054 

2.  05 

0.  036 

10 

1.  70 

0.  068 

2.  05 

0.  026 

11 

2.  20 

0.  062 

2.  05 

0.  048 

12 

2.  10 

0.  086 

2.  00 

0.  060 

13 

2.  20 

0.  096 

2.  00 

0.  094 

14 

2.  30 

0.  110 

2.  00 

0.  092 

15 

1.  90 

0.  116 

2.  00 

0.  092 

16 

1.  85 

0.  138 

2.  00 

0.  110 

10 


TABLE  4 

COMPARISON  OF  OBSERVED  AND  PREDICTED 

T* 

SEARCHER  SPEED:  22  KTS 

TARGET 

SPEED: 

12  KTS 

LAMBDA: 

1. 5/HR 

C(NM.  ) 

D(NM.  ) 

OBSERVED 

PREDICTED 

T* 

P(  T*  ) 

T* 

P(  T*  ) 

20 

3 

0.  70 

0.  036 

0.  70 

0.  036 

4 

0.  60 

0.  056 

0.  70 

0.  046 
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0.  75 

0.  096 

0.  70 

0.  084 
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0.  118 

0.  70 

0.  078 

7 

0.  80 

0.  146 

0.  70 

0.  114 

8 

0.  75 

0.  184 

0.  70 

0.  172 

22 

4 

0.  70 

0.  044 

0.  80 

0.  028 

5 

0.  65 

0.  078 

0.  80 

0.  050 

6 

0.  70 

0.  090 

0.  80 

0.  076 

7 

0.  70 

0.  124 

0.  80 

0.  094 

8 

0.  75 

0.  152 

0.  80 

0.  118 

9 

0.  85 

0.  184 

0.  75 

0.  170 

25 

4 

1.  20 

0.  036 

0.  95 

0.  020 

5 

0.  90 

0.  060 

0.  95 

0.  036 

6 

1.  10 

0.  066 

0.  95 

0.  058 

7 

0.  95 

0.  110 

0.  95 

0.  110 

8 

0.  95 

0.  112 

0.  95 

0.  112 

9 

0.  90 

0.  142 

0.  90 

0.  142 

10 

0.  90 

0.  210 

0.  90 

0.  210 

27 

4 

0.  85 

0.  030 

1.  05 

0.  010 

5 

1.  05 

0.  048 

1.  05 

0.  048 

6 

1.  05 

0.  058 

1.  05 

0.  058 

7 

0.  95 

0.  096 

1.  05 

0.  060 

8 

0.  95 

0.  102 

1.  05 

0.  094 

9 

1.  00 

0.  160 

1.  00 

0.  160 

10 

0.  95 

0.  136 

1.  00 

0.  134 

11 

0.  95 

0.  170 

1.  00 

0.  146 

30 

6 

1.  35 

0.  044 

1.  20 

0.  036 

7 

1.  10 

0.  058 

1.  20 

0.  044 

8 

0.  90 

0.  066 

1.  20 

0.  054 

9 

1.  00 

0.  102 

1.  20 

0.  090 

10 

1.  00 

0.  116 

1.  20 

0.  106 

11 

1.  10 

0.  132 

1.  20 

0.  124 

12 

1.  25 

0.  150 

1.  20 

0.  138 

12 

1.  00 

0.  172 

1.  15 

0.  142 

14 

1.  05 

0.  210 

1.  15 

0.  156 

40 

8 

1.  25 

0.  036 

1.  75 

0.  032 

9 

2.  45 

0.  050 

1.  75 

0.  026 

10 

1.  95 

0.  070 

1.  75 

0.  054 

11 

1.  45 

0.  074 

1.  75 

0.  058 

12 

2.  30 

0.  084 

1.  75 

0.  058 

13 

1.  65 

0.  088 

1.  75 

0.  064 

14 

1.  35 

0.  106 

1.  70 

0.  092 

15 

1.  45 

0.  126 

1.  70 

0.  086 

16 

1.  60 

0.  138 

1.  70 

0.  098 
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